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Abstract 


In PolISAR data processing, deorientation operation is often necessary. The existing deorientation method uni- 
formly deorients all the sub-scatterers of a resolution cell with one orientation angle. For high entropy situation, 
the sub-scatterers have diverse OAs, and the effect of the existing method is unsatisfactory. A novel deorienta- 
tion method is proposed to well treat the high entropy situation. Cloude’s eigen-decomposition to the coherency 
matrix is first carried out. The three eigenvectors are then separately deoriented with their own orientation an- 
gles. Experiments demonstrate that the proposed method is suitable for extraction of urban regions, especially for 


extraction of oriented urban regions. 


1 Introduction 


Deorientation operation plays an important role in fully 
polarimetric synthetic aperture radar (PolSAR) data 
processing. For currently widely used deorientation 
method (the existing method) proposed in [1, 2], the co- 
herency matrix is rotated some angle about the radar 
line of sight (LOS). Only one “mixed” orientation angle 
(OA) is used. However, there are numerous scatterers in 
one resolution cell, and there is a real chance that they 
have different OAs. In oriented urban regions which are 
not aligned in the track direction, there may be roofs 
with azimuth slopes, oriented walls, and oriented dihe- 
drals in one cell [3], for instance. They all have different 
OAs. Only one OA cannot effectively account for this 
situation. A novel deorientation method is thus proposed 
with OA variations (three OAs) taken into considera- 
tion, and the objective of which is to fit the scatterers 
with different OAs in one resolution cell. 


2 Existing Deorientation Method 
2.1 Existing Method 


The existing deorientation method [1, 2] which has 
gained wide applications can be represented as 
1 0 0 
T, = R(@)TR'(6),R(@)=|0 cos20 sin26| (1) 
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where T and T, denote the coherency matrix and the 
deoriented coherency matrix with the existing method 
respectively, and R(@) is the OA rotation matrix. The OA 
0 can be derived by the two methods proposed in [1, 2]. 
The OA is derived as the phase difference between the 


co-polarization terms based on circular polarization co- 
variance matrix in [1], and is equivalently derived by 
minimizing the cross-polarization power in [2]. The de- 


rived OA is 
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where Re(-) denotes the real part of a complex number. 


2.2 Another Interpretation of the Existing 
Method 


Here, the existing method is reconsidered and interpret- 
ed in a different way. Via Cloude’s eigen-decomposition 
[4], the coherency matrix is decomposed as 


3 3 
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where A; k and T; are the eigenvalue, eigenvector, and 
single target, respectively, and 4 = Ay = A, is assumed. 
This provides a statistical interpretation of the target, i.e. 
a three-symbol Bernoulli process is adopted [4]. Thus, 
the existing method of (1) can be rewritten as 
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The three single targets, Tı, T>, and 73, are deoriented 
about the radar LOS with one same OA. Actually as 
mentioned previously, the sub-scatterers in one cell may 
likely have different OAs, and that is to say that Ti, Tə, 
and 7; have their own OAs. Motivated by this point, a 
novel deorientation method is raised in this paper, Le. 
Tı, Tə and T; are separately deoriented with their own 
OAs, respectively. 


3 Proposed Deorientation Method 
3.1 Proposed Method 


The proposed method reconstruct the sub-scatterers of a 
resolution cell, and the sub-scatterers are projected to 
three orthogonal single targets via Cloude-Pottier eigen- 
decomposition of the coherency matrix. Then the three 
single targets are separately deoriented, i.e. 


3 
Tp = ROT R'(G,) (5) 
i=l 
where 7, is the deoriented coherency matrix with the 
proposed deorientation method, 7; is the single target 
obtained with (2), and 6; is the OA of each single target 
T;. 


One noteworthy feature of the proposed method is that 
the real part of T,ı3 element is zero, which is consistent 
with Huynen’s deorientation theory that Huynen param- 
eter H (real part of T,ı3 element) becomes zero after de- 
orientation [5]. But for the existing method, it is not the 
case, the real part of T.23 element is zero. 


The eigenvector k; is here modelled by Kennaugh- 
Huynen’s co-diagonalization of scattering matrix [6, 7], 
and its OA 6, can be afterwards derived. The corre- 
sponding scattering matrix S; of eigenvector k;can thus 
be modelled as 


Si = Rs (9;) Rs (ti )SaiRs (ti) Rs(-G) © 


where Sa is the diagonal matrix whose diagonal ele- 
ments are the complex con-eigenvalues of x; and xp re- 
spectively, and 6; and zq; are the OA and helix angle, re- 
spectively. The unitary transformation matrices are re- 


spectively 
Re(6,) = cosO; -sin 0; (7 
se | sin, cos6; ) 
and 
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The eigenvector k; can be obtained by expanding and 
rewriting the scattering matrix of (6) in Pauli basis as 
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The OA 90; can be derived from (9) as 
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The range of 6; is [-45', 45]. The above derived OA 8; of 
(10) is equivalent to that derived in [8]. 


(10) 


3.2 Physical Interpretation 


The main idea of the proposed method is to rotate the 
three single targets about the LOS separately with dif- 
ferent OAs instead of rotating them together with only 
one OA. We have numerous ways of decomposing one 
distributed target into three single targets, and Cloude’s 
eigen-decomposition is chosen here. As an analytical 
tool, the Cloude’s eigen-decomposition provides a con- 
cise and unique decomposition of the coherency matrix, 
and enables us to interpret the target in a higher dimen- 
sion. Via eigen-decomposition, the coherency matrix is 
decomposed into three statistical single targets (three 
eigenvectors), and then the mean OAs of three single 
targets are all compensated. This results in that the pro- 
posed method is suitable for urban regions extraction, 
especially for characterisation of oriented urban regions. 
In oriented urban regions, there may be roofs with azi- 
muth slopes, oriented walls, and oriented dihedrals and 
so on in one cell. These sub-scatterers are contained in 
the three single targets T, and the proposed method just 
separately rotates the three single targets with their own 
OAs. 


4 Comparisons Between the Two 
Methods 


4.1 Mathematical Comparison 


From the mathematical point of view, the only differ- 
ence between (4) and (5) is that the existing method uni- 
formly rotates the three eigenvectors (the three targets) 
with a same OA 0 about the LOS while the proposed 
method separately and adaptively rotates the three ei- 
genvectors with their own OAs 4), 62, and 6; about the 
LOS. Then the proposed method can be viewed as a 
natural extension of the existing method. The three de- 
oriented eigenvectors by the existing method is still or- 
thogonal, but not the case by the proposed method. The 
two deoriented coherency matrix, i.e. T,and Tp, are both 
positive semi-definite and physical meaningful. 


4.2 Experimental Results 


Here, several experiments are carried out to demonstrate 
the different performances of the two methods. The 
RADARSAT-2 PolSAR data of San Francisco area are 
used [9]. The data are in single-look complex (SLC) 
format. For speckle reduction and getting the coherency 
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matrix format data, a 7X7 refined Lee filter is imple- 
mented on the data. Only a part of the scene is utilized. 


4.2.1 Comparisons Between the OAs 


The OA @ of (2) used in the existing method and the OA 
0, of (10) of the dominant eigenvector used in the 
propposed method are demonstrated and compared, and 
are shown in Figure 1. 
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Figure 1: a) OA of the existing method; (2) OA 6, of the 
proposed method 


(b) 


Four 100X 100 typical patches are selected for further 
comparison. As shown in Figure 2 of Google Earth op- 
tical image, the selected patch 1-4 are ocean region, ur- 
ban region, oriented urban region, and park region, re- 
spectively. 


Figure 2: Google Earth image of selected patches. 


The statistical distributions of the OA 0 used in the 
existing method and the OA 0, used in the propposed 
method for the four patches are shown in Figure 3. 
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Figure 3: Histgrams of the OAs of the four patches. 


From Figure 1 and Figure 3, we can see that 0 and 0, 
are consistent in general, but the OA 0; of the dominant 
eigenvector is noisier as also mentioned in [10]. Figure 
1 shows that the OAs reasonably have large minus val- 
ues in oriented urban region, such as patch 3. 


4.2.2 Impacts of the Two Methods on the 
Model-based Decompsition 


The two deorientation methods have different imapcts 
on the  PolSAR ~~ model-based polarimetric 
decomposition, as will be shown in the folowing. As the 
deorientation methods are the emphasis, original 
Freeman-Durden three-component decomposition 
(FDD) without any modifications [12] is just selected 
here. The FDD is carried out based on the deoriented T, 
and T,. 


First, the negative power problem is tested. It is well 
known that the negative power phenomenon of FDD 
voliates the physical reality. Based on original FDD, 
14.86% pixels have negative power. That numbers for 
FDD with the existing deorientation method and FDD 
with the proposed deorientation method are 8.74% and 
7.77%, respectively. Thus the proposed method is more 
effective in alleviating this physical voliating 
phenomenon. 


Second, the impacts of the two methods on the 
decomposed scattering mechanisms and power 
proportions of the three components are studied. The 
color-coded decomopistion results of FDD, FDD with 
the existing method, and FDD with the proposed 
method are shown in Figure 4, respectively. Blue color 
stands for surface scattering, red color stands for 
double-bounce scattering, and green color stands for 
volume scattering. 


(a) (b) () 


Figure 4: Decomposition results of a) FDD, b) FDD 
with the existing method, and c) FDD with the proposed 
method. 


The average power proportions of the four patches are 
demonstrated in Table 1-4. 


Method Ps Pa 


IE 
EDD | 92.29% | 2.58% | 5.14% 


FDD with the proposed method | 94.17% | 4.43% | 1.4% 
Table 1: Average power proportions of patch 1 ocean 
region. 


FDD with the curret method ee 2.72% | 4.82% 


Method Ps Pa Py 
FDD 29.32% | 51.12% | 19.56% 
FDD with the curret method 31.16% | 52.78% | 16.06% 
FDD with the proposed method | 33.78% | 54.95% | 11.27% 


Table 2: Average power proportions of patch 2 urban 
region. 


Method Ps Pa Py 
FDD -7.718% | 3.47% | 104.3% 
FDD with the curret method 0 27.79% | 72.21% 
FDD with the proposed method | 16.11% | 38.96% | 44.93% 


Table 3: Average power proportions of patch 3 oriented 
urban region. 


Method IR s B d P 


FDD 19.49% | 8.07% | 72.44% 
FDD with the curret method 21.05% | 16.39% | 62.56% 
FDD with the proposed method | 34.72% | 26.70% | 38.58% 


Table 4: Average power proportions of patch 4 park 
region. 


For Figure 4, the red colors in urban regions are 
strenghened by the proposed deorientation method, 
especially in oriented urban regions, such as patch 3. By 
the existing method, the patch 3 oriented urban region is 
not discriminated as double-bounce scattering dominant 
region, but as volume scattering dominant region. By 
the proposed method, this phenomenon is improved, 
and the green colors in patch 3 partially turn to red or 
yellow. The power proportions in Table 2-3 support the 
color change. Thus from the target classification point 
of view, the proposed method is more suitable for urban 
regions extraction, especially for oriented urban regions 
extraction. The green colors in park region, such as 
patch 4, are weakened by the proposed method, but the 
park region can still be correctly discriminated as 
volume scattering dominant region. The power 
proportions in Table 4 support the statement that the 
proposed method gets lower average volume scattering 
power than the curret method. 


5 Discussion 


The proposed deorientation method is designed to better 
treat the high entropy regions. The high entropy regions 
usually include the oriented urban regions and forest 
regions. For oriented urban regions, the proposed meth- 
od has more reasonable and superior performance than 
the existing method as previously demonstrated. For 
forest regions, the proposed method gets lower average 
volume scattering power than the existing method does. 
Although the forest regions can be correctly discrimi- 
nated as volume scattering dominant regions by the 
proposed method, such as patch 4, one question may 
still be raised that does the proposed method underesti- 


mate the volume scattering power in forest regions? 
First, the volume scattering is modelled by the backscat- 
tering from a cloud of dipoles, and the OAs of the di- 
poles obey some probability distribution. The common- 
ly applied uniform distribution and half cosine distribu- 
tion are the centered at zero value, and this is consistent 
with the OA distributions of patch 4 in Figure 3. Se- 
cond, the eigenvector of dipole scattering mechanism 
corresponds to the volume scattering part. Mean OA or 
the center OA of the OA distribution of the eigenvector 
is compensated for each eigenvector by the proposed 
method. As mentioned before, the OAs of forest regions 
usually center at zero generally, therefore, the deorienta- 
tion effect on the eigenvector that corresponds to the 
volume scattering is very small. Thus the proposed 
method does not underestimate the volume scattering 
power for this general case. For the few targets whose 
mean OAs are not centered at zero, the mean OAs of the 
eigenvectors of dipole scattering mechanism will be ro- 
tated, and the volume scattering power from the canopy 
for this case may be underestimated, which may be the 
foreseen limitation of proposed method. However, it is 
difficult to expect that one method can fit all kinds of 
targets. 


6 Conclusion 


Updated from the existing deorientation method 
proposed by Lee et al. and An et al. which utilizes one 
average OA to uniformly rotate the coherency matrix, a 
novel deorientation method whose objective is to fit the 
scatterers with different OAs in one resolution cell is 
proposed. The proposed deorientation method 
separately rotates the three single targets with their own 
OAs obtained by the eigen-decomposion of the 
coherency matrix. One feature of the proposed method 
is that it is consistent with Huynen’s deorientation 
theory, i.e. the real part of the 7,3 becomes zero after 
deorientation, but not the case for the existing method. 
Experimental results demonstrate that the proposed 
method is suitable for urban regions extraction, 
expecially for oriented urban regions extraction. 
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